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FOREWORD 


The  work  reported  in  this  document  was  performed  at  Institute 
of  Seismology,  University  of  Helsinki  unier  the  grant 
OSR-72-2377.  The  author  spent  several  months  at  the  Norwegian 
Seismic  Array  Data  Center  in  Kjeller,  Norway,  and  a month  at 
the  Seismic  Discrimination  Group  in  Lincoln  Laboratory,  M.I.T. 
Cambridge,  benefiting  from  discussions  with  scientists  and 
from  use  of  special  facilities  at  both  places.  Dr.  Clint 
Frazier  of  Lincoln  Laboratory  gave  to  the  author  a useful 
algorithm. 


ABSTRACT 


The  power  spectra  and  waveforms  of  P-wave  sicnals  are  studied 
using  large  aperture  seismic  array  data.  P-wave  pulse  lenghts 
and  attenuation  are  determined  by  deconvolution  of  records 
to  earth  displacement,  and  by  comparing  spectral  shapes  with 
predictions  from  source  models.  A least-square  procedure  is 
used  to  fit  predicted  spectra  to  observations  by  varying  model 
j arameters.  The  attenuation  on  all  ray  paths  from  Asian 
locations  to  NORSAR  studied  was  less  than  t*=0.4,  having  no 
correlation  with  distance  and  a regional  variation.  Larger 
attenuation  is  observed  on  ray  paths  from  SW  North  America  to 
NORSAR.  P-wave  pulse  width  is  defined  as  ratio  between  pulse 
area  and  height.  At  magnitudes  m^  = 4.0,  5.0  and  6.0  the 

estimated  mean  pulse  widths  are  0.4,  0.65  and  1.1  seconds,  &nd 

corresponding  corner  frequencies  are  0.8,  0.5  and  0.3  hz, 

respectively.  The  RMS  variation  corresponds  to  a factor  of 

1.5.  The  explosion  pulse  widths  have  mean  values  0.24  and 
0.31  at  m,  5.0  and  6.0,  respectively,  with  some  disagreement 

D 

in  case  of  NTS  shots.  In  a classification  experiment  using 
statistical  limits  derived  from  NORSAR  data  and  pulse  widths 
measured  at  LASA,  classification  capability  intermediate 
between  the  "spectral"  discriminant  LTM  and  the  n^-Mg 
discriminant  was  achieved,  without  any  regional  corrections. 
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I . INTRODUCTION 


In  this  second  scientific  report  on  study  of  P-waves  observed 
by  the  NORSAR  array,  and  on  seismic  event  classification,  the 
emphasis  is  on  the  determination  of  the  lenght  of  the  P-wave 
pulse  from  small-magnitude  shallow  events.  The  primary  pulse 
from  earthquakes1  and  from  explosions2  is  rather  generally 
accepted  to  be  unidirectional,  save  a possible  small  overshoot 
of  10-20  % which  is  here  usually  disregarded.  The  lenqht  of 
the  pulse  is  an  attractive  quantity  from  the  point  of  view  of 
discrimination,  since  it  is  so  directly  related  to  the  physical 
properties  of  the  source,  and  since  it  should  not  be  so 
strongly  influenced  by  variations  in  wave  attenuation  as 
third  moment  of  frequency  discriminants  and  other  spectral 
discriminants  measuring  energy  at  high  frequencies.  V’e  can 
expect  the  lenght  of  the  wave  pulse  to  be  shorter  for  under- 
ground explosions  than  for  earthquakes. 

For*  example  i-he  source  model  by  Randall"*  predicts  modulation 
of  P- wave  amplitude  spectrum,  caused  by  the  existence  of  two 
acceleration  impulses,  one  at  the  initiation  and  the  other  at 
the  stopping  of  the  P pulse.  In  our  study  it  was  first  tried 
to  observe  spectral  modulation  by  cepstral  methods  . The 
results  were  inconclusive,  and  for  explosions  we  .nave  the 
interfering  pP  echo.  If  we  consider  the  detailed  kinematic 
model  of  e.g.  Savage5  we  find  that  several  stopping  phases  are 
predicted,  corresponding  to  various  phases  of  rupture.  Crustal 
responses  may  also  produce  modulation. 

As  described  in  this  report,  it  is  tried  to  return  to  the 
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am 
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original  earth  displacement  wave  pulse,  since  direct  observation 

of  the  pulse  seems  the  least  ambiguous  method.  This  requires 

correction  for  the  NORSAR  instrument  response.  If  the  pulse 

lenght  is  not  too  large,  it  is  possible  to  get  rather 

accurately  back  to  the  earth  displacement  even  though  who 

sensor  has  no  D.C.  response,  as  shown  in  next  chapter.  It 

seems  also  desirable  to  correct  for  some  estimate  of  wave 

attenuation,  made  practical  by  Carpenter**.  These  corrections 

7 

are  made  using  an  algorithm  by  Frazier  . As  will  be  described, 
it  turns  out  to  be  feasible  to  determine  pulse  lenghts  from 
deconvolved  records. 

Source  duration  has  usually  been  determined  by  spectral 

measurements,  observing  the  high-frequency  cut-off  of  a 

spectrum.  This  is  a more  statistical  method  than  deconvolution 

to  time  domain,  since  by  necessity  it  is  as  much  or  more  a 

measurement  on  the  P coda  as  on  the  primary  pulse.  We  studied 

the  spectra  mainly  to  estimate  wave  attenuation,  which  is  the 

dominating  factor  in  spectrum  at  frequencies  above  2-3  hz.  The 

source  spectral  model  variously  called  as  omega-snuare  (or 

0 

omega-cube)  by  Aki  , or  the  Brune  model,  or  the  spectral 

g 

theory  of  seismic  sources  , is  assumed.  It  is  then  necessary 

to  know  the  corner  frequency  and  the  high-frequency  asymptotic 

slope.  These  were  determined  from  the  data  itself  or  assumed. 

The  remaining  spectral  slope  gives  the  attenuation.  For  the 

explosions,  instead  of  the  "Brune  model",  the  spectral  model 
2 

of  Haskell  with  some  later  modifications  was  used.  It  will 
be  here  referred  to  by  name  "HSB-model" . 

Finally  we  use  our  observations  on  apparent  source  duration 
and  wave  attenuation  in  a classification  experiment  (chapter 
V)  and  for  predicting  values  of  two  other  discriminants 
(chapter  VI) . 
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II. 


MEASUREMENT  OF  PULSE  AREA  AND  WIDTH 


The  teleseismic  P-wave  pulse  is  the  derivative  of  the  source 
time  function  (see  chapter  IV)  and  thus  integration  of  the 
pulse  can  give  information  about  the  source.  In  particular 
the  pulse  area  is  proportional  to  the  source  seismic  moment. 
Since  rhe  area  under  a function  is  equal  to  the  value  of  its 
transform  at  the  origin  and  the  recording  system  has  zero 
response  at  zero  frequency,  objections  can  be  raised  against 
the  possibility  of  determination  of  pulse  area.  However, 
when  the  pulse  width  is  so  narrow  that  at  the  low-freauency 
end  of  the  system  pass-band  it  can  be  approximated  by  an 
impulse,  its  area  can  be  determined.  In  spectral  domain,  it 
is  not  necessary  to  measure  the  power  at  zero  frequency  when 
determing  the  source  seismic  moment  by  spectral  analysis.  A 
measurement  at  any  frequency  situated  within  the  flat  power 
level  surrounding  the  zero  frequency  (as  predicted  by  the 
Brune  theory)  is  sufficient. 

The  greater  the  pulse  width,  the  closer  to  zero  freauency 
one  has  to  go  for  pulse  area  computation.  Thus  it  is  more 
easy  to  measure  areas  for  short  pulses  than  for  Iona  pulses. 

In  the  following  discussion,  by  system  response  is  meant  the 
combined  effect  of  attenuation  and  seismograph  response.  The 
computation  of  these  is,  as  mentioned  in  the  introduction, 
based  on  the  algorithm  by  Frazier  . 

Take  S(f)  to  be  the  spectrum  of  earth  movement.  It  is  assumed 
to  be  the  sum  cf  seismic  noise  E(f)  and  of  a transient  P-wave 
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P(f).  To  study  P(f)  we  want  to  compute  an  approximation 
3e(f)  for  S(f)  from  the  spectrum  T(f)  of  the  recorded  face, 
which  can  written 

T ( f ) = [e  (f ) + P ( f ) ’ R ( f ) + N ( f ) (2.1) 


where  R(f)  is  the  recording-system  response  and  N(f)  is  the 
system  noise.  Dividing  (2.1)  by  R(f)  we  get 

(2.2) 

At  f=0  R (f ) =0 , also  at  very  small  and  large  frequencies 
R ( f ) is  small  and  N(f)/R(f)  can  become  very  larae.  It  is 
necessary  to  limit  estimation  of  S£(f)  to  some  freauency 
band  fL  < t < fH< 


se(f)  = E(f)  + + p(f) 


The  S/N  ratio  decreases  as  the  ranqe  (fL,  fH)  is  widened,  for 

various  reasons.  A possible  way  to  select  values  for  fL  and 

f„  is  to  compute  the  time  transform 
H 


-f. 


se(t)  = 


/ 

-f 


Se(f)e 


iu)  t 


df  + 


H 


/ 

f. 


Se(f)e 


iuit 


df 


(2.3) 


for  various  values  of  fT  and  f„  and  to  select  the  widest 
range  which  gives  acceptable  S/N  ratio  on  the  resulting 
record.  The  limit  f„  is  not  particularly  critical  because 

n 

P(f)  decreases  rapidly  at  high  frequencies.  Cutting  the 
spectrum  at  5 hz  did  have  little  influence  on  our  estimated 
P(t)  but  removed  ringing  problems.  Selection  of  a suitable 
fT  is  critical,  however. 
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We  represent  the  limits  oi  integration  in  (2.3)  by  two 
rectangle  'junctions  n (x),  subtracted  from  each  other 
(figure  II. 1).  The  remainder  function  is 


1(f) 


n (Tf“}  " n (2l^) 


with  transform 


(2.4) 


F ( t ) = 2fp  sir.c(2fHt)  - 2fL  sinc(2fLt)  (2.5) 

where  sinc(x)  = sin  (ttx) /ttx  . F(t)  is  also  shown  in  figure 
II. 1.  For  illustration,  values  fL  = 0.25  hz,  fH  = 5 hz  were 
selected  and  the  time  function  is  shown  as  sampled  with 
sampling  rate  10  samples/sec. 


The  estimated  earth  motion  spectrum  becomes  after  multipli- 
cation with  F (f ) 


se(f)  = H (-jf-)  P(f)  - H (^f-)  P(f)  + F(f)  [E  (f ) 


H 


(2.6) 


N(f) 


It  is  now  assumed  that  either  fH  is  equal  to  the  Nyquist 
frequency  or  that  approximately  P(f)  = 0 at  f > fH-  Then 
the  width  of  the  main  lobe  of  the  function  sinc(2fH0  is 
smaller  than  our  time  resolving  power  or  needs  and  the 
function  can  be  approximated  by  an  impulse  6(t).  We  then 
get  for  the  time  transform  of  3e(f) 


S (t)  = 2f6(t)  * P(t)  - 2 f T sine  ( 2f.  t)  * P(t) 
e H jj  Li 

+ M(t) 


(2.7) 
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Fig.  II. 1.  The  filtering  function  F in  freauency  domain 
(a)  and  in  time  domain  (b) . 


a b 


Fig.  II. 2.  A schematical  P-wave  showing  the  onset 
(a)  and  end  (b)  of  the  P pulse. 


where  M(t)  is  the  transform  of  the  last  (noise)  term  in  (2.6). 


As  current  theories  of  seismic  sources  predict,  the  pulse  is 
assumed  not  to  change  sign.  We  take  it  to  begin  at  t=a  and 
to  end  at  t=b.  We  want  to  determine  the  pulse  area,  ie . the 
integral  of  P(t)  from  a to  b.  Considering  the  convolutions 
in  (2.7)  and  figure  II. 1 we  may  conclude  that  if 


and  noise  and  later  reverberations  can  be  disregarded,  (2.3) 
gives  the  pulse  shape,  apart  from  a base  line  shift.  Berck- 
hemer  and  Burdick  and  Helmberger10  used  this  condition  to 
deconvolve  long  period  seismograms  from  deep  shocks. 


A natural  way  to  measure  areas  of  short  pulses  when  the  base 
line  is  slowly  drifting  is  to  define  the  displacement  at  the 
moment  of  pulse  onset  as  the  base  line  level.  It  is  then 
possible  to  relax  the  condition  fL  <<  (b-a)  1 to  the  less 
demanding  condition  fL  < 0.35 (b-a)-1,  as  will  be  shown  below. 
Noise  M ( t ) will  now  be  disregarded,  since  it  is  stationary 
and  uncorrelated  with  the  signal  its  expected  influence  can 
be  estimated  from  a noise  sample. 


Subtracting  the  amplitude  Se(a)  (the  amplitude  of  Se ( t ) 
immediately  before  the  pulse  onset)  from  both  sides  of  (2.7) 
and  dividing  by  2fH  we  get 


[Se(t)  - Se(a)]  = PU> 
H 


D(t) 


where  the  "distortion  term"  D(t)  is 


(2.8) 
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fL  fI 

D ( t)  = -j—  sinc{2fLt)  * P(t)  - ■£“  sinc(2f^a)  * P(a) 
H H 


or 


fL  b 


D ( t)  - — f P (t  ) [sine  ( 2 f L ( t-T  ))  - sine  (2fL(a-t ) )]  dT 


T rt)  -t 

+ -p  / P(t )[sinc(2fL<t-T) ) - sine (2fL(a-t) )]dx 


(2.9) 


“H  b 


We  investigate  the  first  integral  in  (2.9),  writing  DI(t) 
for  it.  It  is  the  contribution  to  D(t)  by  the  pulse  itself. 
P(t)  is  assumed  to  be  positive  in  (a,b).  If  the  pulse  sign 
would  be  negative,  (2.8)  could  be  multiplied  by  -1. 


Since  t and  t are  in  (a,b)  and 


sine (2fL (a-b) ) ^ sine (2fL ( t-T) ) - 1 
sine ( 2fL (a-b) ) - sinc(2f  ( a—  t ) ) - 1 

one  can  write 

j f r b 

0 < D1(t)  < -^[l  - sine  2f_  (a-b)  ] / P(t)  dx  (2.10) 

fH  L a 

Thus  the  error  caused  by  the  term  DT(t)  in  the  determination 
of  pulse  area  can  be  estimated  from  (2.10)  as 


b f 

0 < / Di(t)  dt  < -j—  (b-a)  [l  - sine  2f_(a-b)] 

a h L 

b 

f P ( T ) dT 

a 


(2.11) 


The  upper  limit  of  the  estimate  (2.11)  for  areal  distortion 
is  plotted  to  figure  II. 3 as  a function  of  the  ratio  between 
the  cutoff  period  T?  = l/fT  and  the  pulse  base  lenght  b-a. 
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pirr  tt  3 The  predicted  maximum  distortion  of  measured 
pi?se  area  as  a function  of  the  ratio  between  cutoff  period 
and  pulse  duration  (the  line)  and  observed  pulse  area 
distortion  (the  circles). 

The  pulse  area  distortion  was  tested  in  the  case  of  an  event 
with  a simple  pulse-like  P-waveform  and  a good  S/N  ratio. 
Deconvolved  records  with  various  values  of  fL  are  shown 
figure  II. 4.  The  value  fH-5  hz  is  used  and  the  attenuation 

amounting  to  t*=0.3  has  been  corrected  for. 


The  event  occurred  at  April  11.  1968,  06-46-27.4  under  the 
Japan  Sea  at  depth  511  km.  Magnitude  mb  is  5.0.  The  width 
b-a  of  the  pulse  base  is  determined  as  1.04  s from  the  low 
f records.  The  per  cent  reduction  in  pulse  area,  relative 
to  the  area  determined  with  fL=0  1 hz,  is  shown  in  figure 
II. 3 as  a function  of  f^  i'-a)'1.  Reduction  of  pulse  area 
follows  closely  the  limit  predicted  by  (2.11).  Oversteppina 
mav  be  caused  by  measurement  inaccuracies.  The  cut-off  at 


T *ACE 


Fig.  II. 4.  A P-wave  record  from  a deep  event  under  Japan  Sea 
deconvolved  using  several  different  long-period  cut-offs, 
shown  at  the  right  margin.  The  original  trace  is  second  from 
top.  Note  the  development  of  a precursor  with  decreasing 
cut-off  period,  and  the  stability  of  shape  at  the  Iona  cut- 
offs. 

[ 

s 
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fT  was  not  quite  abrupt,  taperino  was  used.  However,  the 
discussion  applies  to  any  function  with  transform  roughly  of 
the  sinc(x)  form. 

We  have  not  yet  discussed  the  second  integral  of  (9) . One 
can  give  estimates  for  it,  but  the  distortion  due  to  later 
reverberations  may  or  may  not  be  present,  and  may  be  of 
either  sign.  The  distortion  caused  by  the  pulse  itself  is, 
on  the  contrary,  always  present  and  always  decreases  the 
pulse  area.  It  is  thus  of  more  fundamental  and  important 
nature . 

The  width  of  the  pulse  is  defined  as  the  ratio  area/maximum 
amplitude,  which  is  a standard  definition.  The  maximum 
amplitude  is  also  measured  relative  to  Se(a).  The  implicit 
distortion  on  the  pulse  amplitude  by  DI(t)  term  is  given 
limits  in  (2.10).  Because  both  of  the  pulse  area  and  its 
maximum  amplitude  are  reduced  by  the  term  DI(t),  the  change 
in  their  ratio  should  be  relatively  less  than  the  change  in 
area.  Figure  II. 5 shows  the  apparent  pulse  area  and  width 
as  a function  of  the  cut-off  period/pulse  width  ratio.  The 
area  and  width  are  plotted  relative  to  the  area  and  width 
obtained  using  cut-off  10  seconds.  The  width  is  less  distorted 
by  the  cut-off  than  the  area,  as  suggested  above. 

An  empirical  rule-of-thumb  which  arose  from  deconvolvincr 
a number  of  events  was  that  if  TL/W  (cut-off  period/true 
pulse  width  ratio)  is  greater  than  8,  pulse  area  and  shape 
remain  unchanged,  if  8 < TL/W  < 4 pulse  width  is  unchanged 
but  a precursor  appears,  if  4 < TL/W  < 2 the  pulse  width  is 
reduced  by  the  cut-off,  and  if  TL/W  < 2 no  measurement  is 
possible.  Since  W is  typically  about  half  of  the  pulse 
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Fig.  II. 5 Measured  pulse  width  (triangles)  and  pulse  area 
(circles)  of  the  Japan  Sea  deep  shock,  as  a function  of  the 
cut-off  period/true  pulse  width  ratio.  Both  quantities  are 
shown  relative  to  their  values  measured  with  cut-off  10 
-econds . 


duration  b-a,  these  are  in  agreement  with  the  other  results 
of  this  chapter.  A corollary  of  this  is  that  if  the  measured 
pulse  width  is  less  than  a fourth  of  the  cut-off  period,  we 
can  be  rather  confident  that  no  distortion  has  occured.  If 
the  measured  pulse  width  is  greater  than  that  limit,  apparently 
what  we  can  say  is  that  the  true  pulse  width  is  greater  than 
our  measured  value,  a statement  which  frequently  is  rather 
useful. 

Deconvolution  of  signals  from  explosions  in  Western  Russia, 
at  distances  less  than  30  degrees  from  NORSAR,  gave  usually 
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different  results  than  those  from  more  distant  explosions, 
instead  of  the  simple  waveform  from  distant  explosions  we 
often  observed  complicated  and  elongated  low-frequency 
oscillations  and  in  some  cases  dilatations  of  the  first 
motion.  Our  conclusion  is  that  the  pulse  width  measurement 
by  deconvolution  cannot  be  applied  at  4 < 30°  because  of 
multipa thing  In  the  upper  mantle  (for  multipathing,  see  e.g 
figure  6 of  reference  11)  • 

The  measured  pulse  width  is  influenced  by  the  amount  of 
attenuation  corrected  for  in  the  deconvolution.  This  was 
investigated  by  deconvolving  a few  events  with  several 
different  t'  values.  Results  are  shown  in  figure  II. 6. 


Fia  II. 6.  Pulse  widths  measured  for  explosions  (a)  and 
earthquakes  (b)  after  deconvolution  with  various  values  of 
assumed  t*  Events  are  recorded  at  NORSAR.  Circles  denote 
two  NTS  explosions,  squares  denote  two  E.  Kazakh  shots. 


When  uhe  used  t value  is  increased  from  0.2  to  0.4,  the 
pulse  width  decreases  by  20-30  %.  According  to  these  few 
examples,  for  accurate  determination  of  pulse  width  an 
accuracy  of  0.1  in  t*  estimates  is  sufficient,  and  if  we 
are  willing  to  accept  errors  up  to  50  % the  value  of  t*  has 
no  significance,  as  long  as  it  is  less  than  0.6,  which  we 
believe  to  be  true. 


III.  MEASUREMENT  OF  P-WAVE  ATTENUATION 


Parallerly  with  measurements  in  time  domain,  observations 
in  the  spectral  domain  were  done.  Power  spectral  measurements 
are  of  more  statistical  nature  than  direct  measurements  on 
a pulse.  In  using  the  short  period  sensor  data  for  spectral 
studies  on  P-waves  from  events  in  the  magnitude  ranqe  4.0  - 
6.0  (mb) , one  cannot  avoid  the  influence  of  the  anelastic 
attenuation,  since  the  observed  freauencv  ranne  extends  above 

4 hz.  It  is  tried  to  turn  this  difficulty  into  a benefit  b'f 

+ * 

including  the  attenuation  factor  exp  (-"ft  ),  t as  a free 
parameter,  into  the  spectral  model  which  is  fitted  into  data, 
in  the  hope  of  mapping  attenuation  on  various  wave  paths. 

The  procedure  for  measuring  the  response-corrected  P-wave 
spectra  from  NORSAR  data  was  described  in  the  1st  Scientific 
Report12,  and  will  not  be  repeated  here.  The  average  noise- 
corrected  power  spectrum  of  NORSAR  subarray  beams  is  used. 


We  are  employinq  the  common  sj  ectral  model  discussed  by 
Randall^,  and  proposed  by  Haskell1^,  by  Aki^  and  ty  Brune 
Including  the  attenuation  factor  it  predicts  spectral 
amplitudes 


A ( f ) I = 


A exo(-iTft  ) 
o 


1 ♦ (f/fQ) 


n 


where  A i*  a factor  independent  of  frequency, 
o 

corner  frequency,  usually  associated  with  source 
Both  of  the  factors  exp(-Tft*)  and  (1  + fnfQ  n) 
amplitude  to  decrease  with  increasing  frequency. 


(3.1) 

f is  the 
o 

dimension . 
cav 6 e the 
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If  we  have  available  spectral  information  from  frequencies 
sufficiently  higher  than  the  corner  frequency,  these  two 
factors  can  be  separated  from  each  other.  When  studying 
some  frequency  range,  we  ran  always  determine  the  level  and 
slope  of  the  spectrum.  In  a log  rmpli tude/lir.ear  frecuency 
plot  the  slope  is 


d ( log  1 A 1 ) 
df 


n 

f(f  nf"n+  1) 
o 


TT  t 


(3.2) 


Even  assuming  a value  for  n we  cannot  from  the  slope 

determine  t independently  from  fQ.  But  if  the  frequency 

range  is  sufficient  that  we  can  observe  the  systematic 

2 

change  in  slope,  which  multiplied  by  f is 


f 2 d^ { log | A ] ) M (n-l)n 

_ df5  f n 

O . , 


4- 

then  by  fixing  n we  can  determine  fQ,  and  then  t from 
(3.2)  . 


The  existence  of  modulation  in  P-wave  spectra  makes  it 
necessary  to  observe  a so  wide  frequency  ranae  to  estimate 
the  derivatives  (3.2)  and  (3.3),  that  they  cannot  any  more 
be  described  by  single  average  numbers.  Therefore  the  model 
(3.1)  was  rather  fitted  to  the  logarithm  of  the  observed 
spectra  by  an  iterative  least-squares  procedure.  However 
we  would  note  that  under  conditions  f >>  f (3.2)  and  (3.3) 
become 
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(3.4) 


^ d(loq  A ) _ 

f - at 

f2  d9(log[A|)  _ (n_^n  (3.5) 

When  the  line  described  by  (3.4)  was  fitted  to  the  spectral 
derivative  of  11  earthquakes  in  Central  Asia  in  the  frequency 
range  0.63  - 5.78  hz,  average  values  n = 1.75  +-  0.40  and 
t*  = 0.24  +-  0.03  s came  out.  The  condition  f » fQ  was  not 
well  fulfilled,  however  the  result  for  t*  compares  rather 
well  with  those  described  below. 


Going  now  to  describe  the  results  from  the  iterative  least- 

squares  model  fitting,  the  parameters  determined  in  the 

procedure  were  log(AQ),  log(fQ)  and  t*.  The  parameter  n was 

assigned  values  2 and  3.  Spectra  from  several  events  from  a 

same  region  were  inverted  simultaneously,  keepina  t a common 

parameter  for  all  these  shocks.  To  summarize,  for  each  region 

the  model  was  made  to  fit  observations  from  j events,  each 

having  spectral  estimates  at  i frequencies.  By  inversion  of 

ij  equations,  the  parameters  t , f 0 ( i)  * • • • ' f 0 ( j ) ' Ao(l)""' 

A were  determined. 

o(] ) 

, _ 15 

The  most  probable  values  for  n are  2 and  3.  Recently,  Savage 
suggested  the  value  2.5  for  a circular  fault  with  a rupture 
spreading  at  subsonic  velocity,  and  where  the  slip  at  any 
point  builds  up  gradually.  All  detailed  models  on  body  waves 
generated  by  a deterministic  fault  rupture  predict  oscillations 
in  the  spectrum.  A smooth  spectrum  is  however  useful  as  an 
approximate  modex.  We  are  here  allowing  for  only  one  corner 
£fequgiicy . This  is  to  be  understood  as  an  approximation  for 
the  case  where  the  different  corner  frequencies  are  not  voo 
far  separated.  E.g.  the  slip  time  function  produces  a corner 
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frequency  in  addition  to  the  corner  frequency  associated  with 
rupture  duration.  At  each  corner  frequency  the  exponent  n is 
decreased  by  one^. 


Fig.  III.l.  Estimated  attenuation  as  a function  of  distance, 
when  the  model  parameter  n=2.  Filled  and  open  circles  denote 
ray  paths  from  North  America  and  Asia  t > NORSAR,  respectively 
The  asterisk  denotes  the  path  from  E.  Kazakh  test  site  to 
NORSAR. 

Values  of  t and  of  mean  Q estimated  by  model  fitting  with 
n=2  are  shown  in  figure  III.l.  The  corresponding  results 
with  value  n=3  are  shown  in  figure  III. 2.  Both  results  are 
shown  in  table  UI.l.  In  neither  of  the  figures  the 
attenuation  appears  as  a simple  function  of  distance.  This 
is  not  due  to  inaccuracy  of  measurement,  the  averaae  spectral 
slopes  at  high  frequencies  are  systematically  different 
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between  such  event  groups  as  Sakhalin  and  Kuriles,  or  Greece 
and  Turkey  (see  also  the  plots  in  oh.  Ill  of  ref.  12).  There 
are  only  two  explanations  for  this,  one  is  that  the  earthcuake 
processes  at  the  compared  regions  are  radically  different,  or 
that  the  attenuation  is  different  on  the  respective  ray  paths. 

The  latter  explane  tion  is  much  simpler  and  causes  for  attenuation 
differences  are  not  difficult  to  find  from  regional  tectonics. 


Fig  III  2 Estimated  attenuation  on  ray  paths  as  a function 
* finn  dictance  with  spectral  model  parameter  n 3. 
of  propagation  dista  , act-orisk  denote  ray  paths 

Filled  and  open  c£r®les  and  Eastern  Kazakh  to  NOP.S^R, 

from  western  North  America,  Asia  and  Eastern  k 

respectively.  Explosions  and  earthquakes  were 
different  models. 


A further  check  was  made  on  the  t*  values.  To  all  observed 
earthquake  logo's  were  fitted  a least  square  line,  aivrng 


expected  log^gfo (m^) . For  each  reqion  was  then  separatelv 
computed  the  mean  deviation 

d = AVE[OBS(log1Qfo)  - EXP(log10fo)] 

kt 

Now,  if  the  variations  in  t would  be  caused  by  the 
instability  of  inversion,  such  that  a decrease  in  t is 
offset  by  a decrease  in  f and  vice  versa,  we  would  expect 
to  have  a positive  correlation  between  d and  t*  among  the 
different  regions.  On  the  contrary  there  exists  a negative 
correlation  between  the  two  quantities.  This  can  again 
mean  two  things.  Either  the  corner  frequencies  are 
systematically  higher  at  such  regions  from  where  the  ray 
paths  to  NORSAR  have  a high  Q,  or,  more  probably,  our 
spectral  model  does  not  quite  well  describe  the  reality. 

This  unexpected  correlation  is  observed  both  for  the  cases 
n=2  and  n=3. 

The  RMS  scatter  of  corner  frequencies  around  the  least  squares 

line  is  less  for  the  case  n=3  than  for  the  case  n=2,  by  a 

factor  of  two.  On  the  other  hand,  the  range  of  t*  values  in 

the  case  n=2  agrees  better  with  the  observations  by  Frazier 
17 

and  Filson  and  Noponen  {fig.  III. 13  of  ref.  12).  When 
spectra  rrom  LP  instruments  were  included  in  data,  the  RMS 
values  in  model  fitting  Increased  (events  in  Central  £sia 
and  W.  coast  of  N.  America).  This  is  natural  and  does  not 
mean  that  the  results  would  be  less  reliable,  on  the  contrary 
the  wider  frequency  range  gives  more  reliable  results.  Some 
data  is  shown  with  fitted  models  in  figure  III. 3. 

The  spectral  model  ("HSE-mode1 " ) given  by  von  Seggern  and 
Blandford^ 
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T'ie  explosions  at  E.  Kazakh  have  been  fitted  with  a different  nodel.  The  same  results 
are  written  to  both  of  the  columns  n=2,  n=3. 


(3.6) 


1 

Ao[(2B+l)2  (f-)2  + l]7 


was  fitted  to  the  observed  spectra  from  8 E.  Kazakh  shots, 
using  value  B=2.  The  observed  t*  is  0.20  and  the  resulting 
corner  frequencies  seem  reasonable.  Observed  and  fitted 
spectra  are  shown  in  fig.  III. 4.  The  success  in  fit  was 
apparently  caused  by  the  wide  frequency  range  (up  to  6 hz) 
allowed  by  the  small  attenuation  and  use  of  20  hz  samplinq 
rate.  Analysis  of  other  explosion  regions  where  hiqh 
attenuation  (NTS-NORSAR)  or  our  use  of  data  from  event  tapes 
with  10  hz  sampling  rate  limited  the  frequency  ranqe  to  3 
hz,  gave  corner  frequencies  smaller  than  those  observed  for 
earthquakes,  together  with  relatively  low  t*  values.  The 
reason  is  probably  the  strong  modulation  of  spectrum  by  the 
pP  phase,  requiring  a wide  frequency  range  in  the  estimation 
of  the  general  trend  of  spectrum.  Direct  comparison  between 
NTS  and  E.  Kazakh  explosion  spectra  recorded  at  NORSAR  leaves 
no  doubt  of  the  much  stronger  attenuation  on  the  NTS-NORSAR 
path,  the  suggested  difference  being  about  0.28  in  t* 

( figure  II 1 . 5) . 

If  we  add  0.28  to  the  absolute  t value  determined  for  the 
E.  Kazakh-NORSAR  path,  we  get  0.48  as  the  absolute  t*  on 
NTS-NORSAR  path,  in  a reasonable  aareement  with  the  value 
0.41  obtained  by  Filson  and  Frazier  . 

A further  experiment  was  to  compute  the  averaae  spectrum  of 
4 NTS  shots  with  mean  111^=5.2,  and  the  average  spectrum  of 
6 E.  Kazakh  shots  with  mean  1^=5. 5,  to  restrict  the  corner 
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Fia  HI  3 Observed  spectra  from  a few  earthquakes  in 
Central* Asia,  together  Sith  spectra  predicted  by  the  Brun 
model,  fitted  into  the  observations  by  varying  parameters 
several  least-squares  iterations.  The  crosses  at  the 
indicate  the  position  of  fQ  and  the  standard  deviations  o 
f and  A , for  each  curve. 


Fig.  III. 4.  Observed  spectra  from  E.  Kazakh  explosions, 
together  with  spectra  predicted  by  the  Haskell  - von  Seqqern  - 
Blandford  model  fitted  into  observations  by  varyinq  parameters 
in  successive  least-squares  iterations.  The  crosses  close 
to  1 hz  indicate  the  position  of  f and  the  standard  deviations 
of  fQ  and  Aq,  for  each  model  curve? 
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frequencies  of  the  two  average  spectra  as  equal,  and  to 
determine  the  one  f and  the  two  t*  in  a simultaneous 
inversion  to  the  HSB-model.  Results  are  given  in  table 
III  .2 . 


Fig.  III. 5.  Spectral  ratio  between  NTS  and  E.  Kazakh 
explosions  in  the  magnitude  range  (nv)  5. 0-6.0.  Assumption 
of  similar  average  source  spectrum  gives  for  attenuation 
difference  between  paths  NTS-NORSAR  and  E.  Kazakh-NORSAR  the 
value  t*=0.28,  from  the  visually  fitted  line  in  the  fiaure. 
This  figure  is  in  agreement  with  the  attenuation  difference 
between  southwestern  North  America  - NORSAR  and  Central  Asia 
NORSAR  paths,  determined  from  the  ratio  of  earthquake 
spec tral2  as  0.279  i 0.004. 


Table  III. 2 results  of  an  inversion  experiment. 

t* (NTS-NORSAR)  0.403  - 0.076 

t*(E.  Kazakh-NORSAR)  0.150  - 0.027 

The  common  f 0.61-0. 25  hz 


DUO 


2.5  dB 


and  the  observed  and  predicted  spectra  are  shown  in  figure 
III. 6.  In  averaging  spectra  from  different  shots  the  effect 
of  the  pP  modulation  probably  decreased,  and  the  high- 
frequency  E.  Kazakh  spectra  stabilized  the  common  solution. 
The  impression  from  the  figure  is  that  the  lowest  freauency 
band  in  the  NTS  spectral  average  has  a disproportionally 
high  amplitude,  distorting  the  solution  so  that  the  common 
f approaches  the  earthquake  fQ  at  a similar  magnitude. 


Fig.  III. 6.  Average  spectra  for  4 NTS  explosions  (filled 
circles)  and  6 E.  Kazakh  explosions  (open  circles)  and  HSB 
model  spectra  simultaneously  fitted  to  both  data.  The  dB 
scale  is  relative  to  the  zero-f reauency  level  of  the 
theoretical  curves. 
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IV.  DURATION  OF  EARTHQUAKE  AND  EXPLOSION  SOURCES 


At  this  stage  no  interpretation  of  the  observations  in  terms 

of  a finite  source  model  will  be  done.  Only  the  apparent 

source  time  function  S(t)  will  be  discussed.  If  a point 

source  consisting  of  force  couples  is  assumed,  S(t)  becomes 

the  time  history  of  the  forces,  the  radiated  teleseismic  P- 

wave  being  the  derivative  dS(t)/dt.  The  difference  between 

earthquakes  and  explosions  lies  then  only  in  the  geometry  of 

the  force  couples,  a double  force  couple  describing  the 

earthquake  source19,  three  perpendicular  couples  without 

20 

moment  the  explosion  source 


A point  source  is  in  principle  an  unsuitable  model,  since  the 
wavelenghts  of  the  observed  waves  are  of  the  order  of 
earthquake  size.  However,  since  our  observations  for  each 
shock  are  limited  to  a single  direction,  there  will  be  no 
need  to  use  finite  source  models  giving  different  S(t)  to 
different  directions. 

As  S ( t)  is  monotonous1'21,  save  for  a possible  small  over- 
shoot2 which  we  disregard  here,  the  observations  in  time 
domain  and  frequency  domain  can  be  conveniently  related 
through  the  reciprocal  concepts  of  spectral  width  and  pulse 
width  . Defining  the  width  of  a "unction  as 


we  get  for  the  width  of  the  amplitude  spectrum  of  the  Brune 
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model  the  value  nf  , and  for  the  corresponding  model  with 

° -3 

the  high-frequency  spectrum  proportional  to  f the  value 

3_1/2nfo.  If  the  pulse  shape  in  time  domain  is  an  even 

function  the  pulse  width  is  the  inverse  of  the  spectral 

width.  Thus  the  relationships  between  the  pulse  width  Wp 

and  the  corner  frequency  are 


(4.2) 

(4.3) 


for  the  values  n=2  and  n=3,  respectively. 


Pulse  width  measurements,  made  by  the  methods  described  in 
chapters  II  and  V,  are  plotted  in  figure  IV. 1 against  the 
corner  frequencies  measured  as  described  in  chapter  III,  for 
the  same  events.  Results  from  both  of  the  assumptions  n=2 
and  n=3  are  shown.  Only  events  with  distance  from  NORSAR 
greater  than  30  degrees  are  included. 

When  comparing  the  lines  predicted  by  (4.2)  and  (4.3)  to 
the  plotted  data,  one  observes  than  in  case  of  n=2  about  2/3 
of  the  points  fall  reasonably  close  to  the  prediction.  Among 
the  1/3  which  have  fQ  over  twice  the  value  predicted  from 
the  pulse  width  we  find  all  of  the  7 plotted  Kuriles  shocks, 

4 of  the  19  Sakhalin  shocks  and  3 of  the  14  Central  Asia 
shocks.  In  case  of  n=3  such  ipurious  high  frequencies  do 
appear  only  for  a few  of  shocks. 

At  least  a couple  of  explanations  suggest  themselves.  In 
chapter  III  it  was  noticed  that  for  the  ray  paths  with  high 
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^iq  IV  1.  Measured  corner  frequencies  plotted  against  the 
neasured  pulse  widths  for  the  two  cases  n=2  and  n=3  The 
:ontinuous  lines  show  the  expected  relations.  Filled  circles 
ienote  E.  Kazakh  explosions. 


ft 

t values  also  the  corner  frequencies  were  above  average.  It 
was  then  suggested  that  either  the  attenuation  model  or  the 
source  spectral  model  is  not  accurately  describing  the  reality, 

ft 

causing  a coupling  between  the  parameters  t and  f . This 

could  explain  why  all  the  Kuriles  shocks  have  excessive  f 

- 1 ° 

when  n=2.  We  can  also  have  a,  say,  f asymptotic  spectral 

dependence  up  to  a fairly  high  second  corner  frequency  in 

case  of  some  events.  This  seems  to  apply  to  the  event  "229" 

in  Central  Asia,  for  which  the  spectral  slope  turns  normal 

around  4 hz . The  pulse  lenght  needs  not  then  to  be  particularly 

short,  but  the  pulse  has  "sharp  corners".  In  case  of  the 

event  "105"  we  have  after  i_he  primary  pulse  (with  duration 

1.4  s and  width  1.1  s)  a large-amplitude  oscillation  with 

frequency  1 hz,  which  probably  produces  the  sharp  corner 

observed  in  spectrum  at  1 hz  (spectrum  shown  in  reference  12, 

fig.  II. 1).  Such  la.er  reverberations  can  be  caused  by  a 

22 

thick  low- velocity  layer  above  a shallow  source 

We  make  a preliminary  assumption  that  the  high  corner 
frequencies  which  disagree  with  corresponding  pulse  widths 
are  a result  of  our  crude  source  spectral  model  not  always 
being  able  to  model  the  spectrum  correctly,  and  that  they 
need  not  imply  very  short  widths  for  the  primary  pulse. 

In  figure  IV. 2 the  90  % confidence  limits  for  the  t: ue  pulse 
23 

width  from  various  sets  of  data  are  shown  as  functions  of 
m^.  The  following  assumptions  have  been  made:  the  width/ 

magnitude  relationship  is  roughly  linear  in  the  magnitude 
interval  of  study,  the  individual  width  measurements  are 
distributed  randomly  and  roughly  normally  around  the  parent 
width/m^  line,  the  measured  corner  frezuencies  can  be  trans- 
formed to  pulse  widths  by  (4.2)  and  (4.3),  and  the  errors  in 
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magnitude  are 
width. 


much  smaller  than  the  scatter  in  pulse 


pl„  IV  2 The  90  % confidence  limits  for  true  pulse  width 

r Us  s« 

explosions  (the  lower  ^urY®®dmf  m corner  frequency  measurements/ 
cfherliirtf arff^m  direct  measurements  on  pulse  width. 


The  different  samples  which  are  described  by  fig.  IV. 2 are: 

(1)  Corner  frequencies  measured  under  assumption  n=3  from 

spectra  of  91  earthquakes  recorded  at  NORSAR,  and  transformed 

to  pulse  widths  by  (4.3).  The  it^  range  is  4.1  to  6.0,  limits 

are  plotted  by  dotted  line.  (2)  Pulse  widths  measured  from 

deconvolved  NORSAR  records  of  53  earthquakes,  forming  a subset 

of  previous  the  set.  Magnitude  range  is  4.1  to  6.0,  limits 

are  plotted  by  broken  line.  (3)  Pulse  widths  measured  from 

deconvolved  LASA  records  of  59  earthquakes24,  in  magnitude 

range  4.1  - 6.1.  Limits  are  plotted  by  continuous  line. 

(4)  Corner  frequencies  measured  using  the  HSB-model  from 

spectra  of  8 E.  Kazakh  explosions  recorded  at  NORSAR.  The 

magnitude  range  is  4.2  - 5.8,  limits  are  plotted  with  dotted 

line.  (5)  Pulse  widths  measured  from  deconvolved  NORSAR 

records  of  22  explosions  in  Kazakh  and  NTS.  Magnitude  range 

is  4.4  to  5.8,  limits  are  plotted  with  broken  line.  (6)  Pulse 

widths  measured  from  deconvolved  LASA  records  of  10  explosions 

24 

at  various  locations  , the  magnitude  range  is  4.8  to  6.1. 
Limits  are  plotted  with  continuous  lines. 

Figure  IV. 3 seems  to  give  a consistent  picture  of  the  P-wave 
pulse  width  and  its  change  with  magnitude  for  earthquakes. 

For  the  explosions  there  is  a disagreement  between  the  data 
sets  (5)  and  (6),  at  magnitudes  below  5.5.  It  does  not  seem 
difficult  to  clarify  the  situation  with  further  research. 

Of  more  interest  from  the  point  of  view  of  study  of  earthquake 
mechanism  would  be  the  plot  of  pulse  width  against  source 
seismic  moment,  which  can  be  computed  from  the  pulse  area. 

Such  a plot  will  however  be  worse  as  an  identification 
diagram,  due  to  the  small  moments  of  explosions.  The  seismic 
moments  could  also  be  plotted  against  the  magnitude  m^,  giving 
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a picture  resembling  the  m^/Mg  P^ot»  since  at  - 6 Mg  is 
proportional  to  the  moment25.  The  influence  of  source  depth 
would  be  removed,  however,  and  the  influence  of  source 
mechanism  would  be  different.  A preliminary  plot  snowed  that 
a moment/n^  plot  gave  no  better  discrimination  than  a width/ 

m^  plot. 

The  union  and  intersection  of  the  confidence  regions  of 
figure  IV. 2 are  shown  in  figure  IV. 3 by  the  hatched  and  black 
regions,  respectively.  Earthquakes  and  explosions  are  treated 
separately.  The  confidence  limits  computed  from  the  fQ 
derived  under  assumption  n=2  are  also  shown,  by  broken  lines. 
The  region  bounded  by  these  limits  has  an  intermediate  position 
between  the  explosion  and  earthquake  90  % confidence  regions 
derived  from  the  other  sets  of  data.  As  mentioned  above,  our 
assumption  was  that  spurious  high  corner  frequencies  have 
distorted  the  average  properties  of  this  set  of  data. 

Comments  on  NTS  explosions.  Even  after  taking  account  the 
large  attenuation  on  the  ray  path  NTS-NORSAR,  the  pulse 
widths  of  NTS  explosions  seem  systematically  somewhat  larger 
than  those  of  the  Eastern  and  Central  Kazakh  explosions,  e.g. 
causing  the  disagreement  between  the  data  sets  (5)  and  (6) . 

One  explanation  for  this  could  be  that  on  the  frequency  1.3 
hz,  around  which  the  mb  is  measured,  an  attenuation  difference 
of  .28  in  t*  produces  a factor  of  3 difference  in  spectral 
amplitude.  This  rough  calculation  suggests  that  for  two 
explosions  of  similar  yield  fired  under  similar  conditions 
at  the  two  sited,  the  NTS  explosion  would  be  given  up  to  0.5 
units  smaller  magnitude  than  for  the  E.  Kazakh  explosion,  at 
NORSAR.  If  the  attenuation  conditions  on  ray  paths  from  tne 
two  sites  to  the  majority  of  the  stations  are  similarly 
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different  than  on  paths  to  NORSAR,  the  NTS  explosion  magnitudes 
should  be  corrected  upwards  by  some  source  factor,  for  comparison 
The  disagreement  between  Kazakh  and  NTS  pulse  widths  would 
then  decrease. 


Fig.  IV. 3.  The  union  (hatched)  and  intersection  (black) 
of  the  confidence  regions  shown  in  previous  figure,  separately 
for  earthquakes  and  explosions.  The  confidence  limits  for 
an  earthquake  set  analysed  in  spectral  domain  under  assumption 
n=2  are  shown  by  broken  lines.  Our  interpretation  is  that 
they  are  distorted  by  spurious  high  corner  frequencies. 


Figv."  xll.6  suggests  excessive  low  frequency  energy  for 
E Kazakh  explosions.  In  deconvolving  the  NTS  records  from 
NORSAR,  sometimes  slow  motions  of  large  amplitude  were 
observed.  We  cannot  give  any  explanation  for  them,  we  are 
not  even  sure  of  their  reality.  But  inevitably  they  distort 
our  spectral  measurements.  In  figure  IV. 4 are  shown  the 
NORSAR  beam  record  and  deconvolutions  with  cut-offs  10  and 
2.3  seconds  for  the  NTS  explosion  on  day  156  of  1973,  time 
17-11-32.5,  displaying  a large  later  pulse  with  duration  of 

appr . 4 seconds . 


f.  0.45 


Fig.  IV.  4.  NORSAR  beam  trace  (top)  of  a ,NTS . ®xpl°s^' 21  ? 
deconvolution  with  cut-off  10  seconds  and  with  cut-off  2.3 

seconds  (botto'.a)  . 


V.  CLASSIFICATION  EXPERIMENT  USING  PULSE  WIDTH 


The  pulse  width  determination  method  developed  using  the 
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NORSAR  data  was  tried  on  a set  of  data  collected  by  Lacoss 
from  LASA  records  on  1966-1968.  It  was  possible  for  the 
author  to  use  the  console  system  of  the  Seismic  discrimination 
group  in  Lincoln  Laboratory,  M.I.T.  The  system  was  well 
suited  for  this  kind  of  analysis  and  speeded  up  the  work. 

The  procedure  followed  is  now  described: 

(1)  The  LASA  beam  record  was  displayed  on  the  screen  and  by 
varying  the  gain,  the  onset  time  and  polarity  were  observed. 

(2)  The  beam  was  deconvolved  with  t*=0.3  s and  using  low- 
frequency  cut-off  periods  of  10,  7,  4,  2.5  and  1.75  sec. 

The  resulting  five  traces  were  displayed  on  the  screen. 

(3)  Generally  the  S/N  ratio  decreases  with  increasing 
cut-off  period.  The  trace  with  longest  cut-off  period  and 
with  acceptable  S/N  ratio  was  selected.  It  was  checked  that 
the  pulse  had  a correct  polarity  and  began  appr.  0.3  seconds 
prior  to  the  onset  on  beam.  If  no  sufficient  S/N  ratio  was 
observed  or  there  was  a disagreement  in  pulse  onset  time  or 
polarity  with  the  beam  record,  the  event  was  usually  placed 
to  the  category  of  non-measurable  events. 

(4)  The  pulse  width  was  measured  semi-automatically . 

Two  positions  were  manually  defined  on  the  selected  trace. 

The  first,  t^,  was  the  expected  onset  time  on  deconvolved 
record,  the  second,  t2»  was  approximately  on  the  pulse 
maximum.  Then  an  integration  procedure  was  triggered,  giving 
the  pulse  width  as 
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where  t is  the  point  where  min  A(t)|t^  ^ t < t2  is  observed, 
t is  min  ( t , t.  ),  where  t is  the  first  point  where  A(t)  < A(t  ) 

e 3 D 3 

and  t > t2,  snd  tfa  is  the  first  point  where  dA(t)/dt  > 0, 
t > t2  snd  A ( t)  <0.3 (max  A (t) J tg  < t < t) . Am  = max  A(t) | 

fcs  < < fce  * 

In  figure  V.l  deconvolved  records  of  a magnitude  mb=5.2 
earthquake  in  western  Turkey  are  shown.  Pulse  width  was 
measured  from  the  record  with  cut-off  period  7 seconds.  The 
long  and  short  vertical  lines  show  the  manually  adjusted  times 
tl  and  t2.  The  pulse,  limited  by  times  tg  and  te,  is  shown 
as  the  top  trace.  The  computed  pulse  width  is  1.28  seconds, 
giving  a cut-off /width  raeio  of  5.5  which  according  to  the 
discussion  in  chapter  II  indicates  that  pulse  width  was  not 
reduced  by  cut-off.  The  duration  of  the  pulse  is  2.24 
seconds.  The  shown  event  illustrates  an  interesting  pheno- 
menon, the  pulse  has  two  maxima.  Incidentally,  when  the  cut- 
off period  decreases  the  two  maxima  separate  to  two  pulses. 

A pulse  width  measurement  from  the  trace  with  cut-off  2.5 
would  give  too  small  pulse  width,  in  itself  a natural  result 
from  too  small  true  width/cut— of f ratio,  but  in  this  case 
unobservable  since  the  minimum  in  the  pulse  disrupts  the 
rule-of-thumb  given  in  chapter  II.  Even  the  underestimated 
pulse  width  (0.45  s)  falls  within  the  earthquake  population 
in  the  classification  diagram  V.2,  however. 


Pulse  widths  measured  for  the  first  84  events  of  the  LASA 

A Jl 

data  tape  created  by  Lacoss  are  shown  in  figure  V.2.  The 
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rule  that  a pulse  width  is  measured  first  and  the  nature  of 
the  event  checked  only  afterwards  was  strictly  obeyed.  In 
addition  a few  confidence  limits  and  the  earthquake  and 
explosions  mean  square  pulse  width  lines,  computed  from 
independent  samples  of  pulse  widths  measured  from  NORSAR 
records,  are  drawn  to  fig.  V.2. 


Fig.  V. 1 . Deconvolved  records  from  A mj=5.2  earthquake  in 
western  Turkey,  recorded  by  LASA.  The  figures  at  right 
margin  give  the  cut-off  period  for  each  trace.  The  top  trace 
shows  the  semi-automatical ly  selected  pulse  from  the  record 
with  cut-off  7.0.  The  times  t ^ and  t2  were  manually  set. 
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Fig.  V.  2 . Pulse  widths  measured  for  earthquakes  (open 
circles)  and  explosions  (filled  circles)  recorded  by  LASA. 
The  statistics  plotted  to  the  figure  are  derived  from 
independent  earthquake  and  explosion  populations  recorded 
by  the  NORSAR  array.  LSQ=least  squares  line  for  earthquakes 
95Q=95%  confidence  limit:  for  a single  earthquake  pulse 

width,  R1Q=99%  confidence  limit  for  detecting  outliers  from 
the  earthquake  population.  LSE,  95E,  R1E  corresponding 
quantities  for  explosions.  DL=preliminary  decision  line. 


The  sample  of  earthquakes  used  in  computing  the  statistics 
is  the  set  (5)  of  chapter  VI,  it  consists  of  53  quakes,  ^he 
explosion  sample  consists  of  17  Eastern  and  Central  Kazakh 
shots,  the  NTS  shots  were  excluded  in  computing  the  statistics, 
since  they  seemed  to  belong  to  a different  population.  No 
NTS  shots  are  included  in  the  LASA  data  tape.  It  is  necessary 
later  to  unify  all  observed  explosions  under  a single 
description,  to  be  more  confident  that  no  anomalous  explosion 
populations  will  appear,  to  which  our  expectations  do  not 
hold. 
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The  computed  statistics  are  : (1)  The  least  squares  lines 

for  NORSAR- recorded  earthquakes  (LSQ)  and  explosions  (LSE) . 

(2)  The  one-sided  95  % confidence  limits  for  a single  sample, 
lower  limit  for  earthquakes  (95Q),  upper  limit  for  explosions 
(95E) . (3)  The  99  % confidence  limits  for  rejecting  outliers 

from  the  explosion  population  ( R1E)  and  from  earthquake 
population  (R1Q) . A preliminary  decision  line  (DL)  is  also 
drawn  to  the  figure.  It  divides  the  distance  between  the 
least-square  lines  in  the  ratio  of  the  residual  RMS  values 
of  explosions  and  earthquakes 

The  small  scatter  of  the  original  Kazakh  explosion  pulse 
width  data  at  NORSAR  is  reflected  in  the  outlier  test  for 
explosions,  which  suggests  all  except  3 of  the  quakes  not 
to  belong  to  the  explosion  population,  with  a 99  % confidence. 
On  the  contrary  the  large  scatter  of  earthquake  pulse  widths 
causes  that  none  of  the  explosions  when  treated  singly  is 
suggested  to  be  an  outlier  from  the  earthquake  population 
at  the  same  level  of  confidence.  It  is  thus  easier  to 
identify  quakes  than  explosions. 
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One  of  the  earthquakes  falls  out  of  the  95  % confidence 
region  for  single  earthquakes,  and  into  the  95  % confidence 
region  for  explosions.  That  is  a deep  shock  in  Rumania 
(depth  151  km) . There  are  other  deep  shocks  with  normal  pulse 
widths,  however.  Two  other  earthquakes  are  marginal,  they 
are  shallow  shocks.  Results  are  summarized  in  table  V.l. 

Table  V.l.  Results  of  the  discrimination  experiment 


source  type 

classified  as: 
explosion  earthquake 

marginal 

rejected 

total 

EXPLOSION 

10 

0 

0 

1 

11 

EARTHQUAKE 

1 

56 

2 

14 

73 

The  one  rejected  explosion  was  clipping.  In  figure  V.3  are 
shown  the  deconvolved  records  of  9 explosions,  and  of  the 
three  earthquakes  not  classified  as  earthquakes.  We  think 
that  all  of  these  quakes  are  in  some  way  different  from  the 
explosion  records.  Q/l  is  the  deep  shock,  it  has  no  pP 
reflection  following  the  primary  pulse,  as  all  the  explosions 
have.  Q/2  is  the  most  problematic  of  all.  It  is  a shallow 
shock  with  magnitude  1^=4. 7 at  Kuriles,  under  the  ocean.  A 
secondary  negative  impulse  appears,  resembling  the  pP 
reflection  from  explosions.  The  later  reverberations  of  the 
signal  are  however  larger  and  with  a longer  duration  than  for 
any  of  the  explosions.  The  S/N  ratio  is  bad.  The  event 
Q/3  is  a nifa=6  shock  in  SE  Europe.  Apparently  we  are  measuring 
a small  precursor  of  a large  and  complex  shock.  The  appearance 
of  the  record  is  clearly  not  explosion-like. 

The  original  beam  traces  of  the  same  12  events  are  shown 
in  figure  V.4.  The  amplitude  gain  of  each  record  is  random. 
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Fig.  V.3.  Deconvolved  records  for  the  explosions  and  the 
three  earthquakes  not  classified  as  earthquakes. 

E=Eastern  Kazakh,  W=Western  Kazakh,  A=Aleutians.  Q stands 
for  quake.  Code  numbers  refer  to  appendix.  Amplitudes 
are  normalized. 
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The  polarity  of  the  event  Q/2  is  somewhat  uncertain,  during 
analysis  it  was  decided  to  be  positive.  A supplement  to  the 
pulse  width  discriminant  is  the  polarity,  since  the  determination 
of  the  width  requires  the  pulse  polarity  to  be  first  determined. 
At  distances  above  30  degrees  all  explosion'.,  from  which  the 
pulse  width  can  be  determined  display  unmistakably  compressions. 
For  earthquakes  also  dilatations  are  observed. 

Nexc  the  95  % one-sided  confidence  limits  for  a single  sample, 
the  least  square  lines  and  th°  preliminary  decision  line  are 
computed  from  the  LASA  event  sample  itself.  These  are  thus 
statistics  which  can  be  used  in  classifying  new  LASA- recorded 
events.  They  are  shown  in  figure  V.5.  The  data  points  are 
also  plotted  to  the  figure. 

In  table  V.2  are  shown  the  number  of  events  accepted  to  the 
"measurable"  category  and  th®  total  number  of  analysed  events, 
and  the  percentage  of  accepted  measurements  in  several 
magnitude  categories. 

Table  V.2.  Statistics  on  measurement  capability 


m^  class 

total 

accepted 

acc.  % 

4. 0-4. 9 

33 

2 L 

63% 

ON 

• 

m 

l 

o 

• 

in 

43 

f,2 

98% 

over  6.0 

7 

61 

1)  one  of  the 

large  explosions 

clips . 

The  pulse  width  measurements  can  often  be  extended  to  such 
explosions,  for  which  Mg  is  not  measurable  and  m^-Mg  classifier 
cannot  be  used. 
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5 
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Fia  V 5.  Statistics  computed  from  the  LASA  pulse  wioth 
da?;  sample  analysed  in  this  chapter.  Symbols  are  the 

same  as  in  figure  V.2.  The  widtn  ”ea^e^n^ea^s^^n 
plotted.  These  statistic*  were  nor  used  in  the  discri 

mination  experiment. 
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VI.  IMPLICATIONS  TO  OTHER  DISCRIMINANTS 


The  knowledge  on  source  duration  and  P-wave  attenuation  which 
has  been  achieved  in  the  previous  chapter  may  be  applied  to 
exptrapolate  and  predict  values  of  such  other  discriminants 
which  are  influenced  by  these  quantities. 

The  LTM  discriminant  measures  the  amount  of  P-wave  high 
frequency  energy  relative  to  the  low  frequency  energy.  The 
normalized  third  moment  of  the  noise-corrected  power  spectrum 
around  the  frequency  origin  is  computed.  By  defining  the 
quantities  t*  and  fQ,  and  applying  Brune  and  the  HSB  model 
for  earthquakes  and  explosions,  respectively,  sets  of  predicted 
spectra  were  computed.  These  were  then  multiplied  by  the 
NORSAR  amplitude  response  and  LTM  values  were  computed. 

In  figure  VI.  1 the  results  are  compared  with  earthquakes  and 
explosions  in  Asia  and  North  Am  rica.  t values  from  figure 
III.l,  earthquake  fQ  values  derived  from  figure  IV.  2 and 
explosion  fQ  values  derived  from  figure  V.2  using  (4.2)  were 
used.  Agreement  is  acceptable  for  earthquakes  and  explosions 
in  Central  Asia,  Turkey,  and  Sakhalin.  The  predicted  curves 
lie  too  high  for  North  America  (t  =0.5)  and  too  low  for  Kuriles 
(t*=0.2).  It  is  again  observed  that  even  using  our  estimated 
attenuation  figures,  events  in  North  America  have  still 
weaker  high  frequency  energy  than  expected,  and  events  in 
Kuriles  have  still  stronger. 

As  mentioned  in  chapter  III  the  probable  cause  for  these 
phenomena  is  a spectral  model  which  does  not  quite  well  fit 
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t reality.  One  possible  consequence  is  that  attenuation  on 
North  Air.erica-NORSAR  paths  is  underestimated  and  on  Kuriles- 
NORSAR  path  overestimated.  Another  observation  is  that  the 
observed  earthquake  LTM  values  decrease  more  steeply  with 
increasing  magnitude  than  the  predicted  values. 


Fig.  VI.  1.  At  NORSAR  observed  LTM  values  from  ref.  12  compared 
with  predicted  values.  Explosion  and  earthquake  observations 
are  plotted  with  filled  and  open  symbols,  respectively. 
Continuous  and  broken  lines  denote  predictions  for  explosions 
and  quakes.  (a)  Events  in  south-western  North  America. 

(b)  Explosions  in  western  Russia,  earthquakes  in  southern 
Greece (circles)  and  Turkey  (triangles).  (c)  Events  in  Central 
Asia.  (d)  Earthquakes  at  Kuriles  (circles)  and  Sakhalin 
(triangles) . Almost  all  of  the  events  in  the  figure  are 
shallow,  and  none  is  deeper  than  IOC  km. 
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Next  it  will  be  attempted  to  combine  our  results  with 
mb“Ms  observations • The  magnitude  Mg  is  first  derived. 

The  amplitude  of  Rayleigh  waves  at  period  20  seconds  is 
given  as 

AR20  7?  Po  DR20 WR(h» S,E)  | 

for  a step  source  time  function.  PQ  is  seismic  source 
moment,  dr2q(A)  Is  t^le  luence  of  distance  on  amplitude, 

W is  the  effect  of  source  depth  h,  mechanism  S and  Earth 

**  26 
structure  E.  Douglas,  Hudson  and  Kembhavi  showed  that  a 

separation  of  factors  in  this  way  is  possible.  Since  M 

s 

is  defined 

Ms  * lo9  <W  + V4) 


where 


Q ( A)  + log  (Dr2Q(A)  ) = eR 

en  being  a constant,  related  to  the  origin  of  the  M scale 

K S 

w’  get  for  M 

s 

Ms  = log  7?  + log  p0  + eR  + log |wR(d,S,E) | 

For  the  magnitude  one  can  give  t'e  following  derivation 
The  amplitude  spectrum  of  a P wave  is 

Ap(f)  = F ( f ) PQ  Dp  ( A,  h,  E ) exp(-TTft*)  |wp(S,E)  | 

where  F(f)  is  the  amplitude  spectrum  of  the  source  time 
function,  Dp  stands  for  geometrical  attenuation  and  the 
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effect  of  the  free  surface  at  receiver  end,  t is  the  ray 
path  attenuation  and  to p gives  the  effect  of  the  properties 
of  the  supposed  homogeneous  medium  surrounding  the  source, 
and  of  the  source  radiation  pattern,  assumed  to  influence 
only  the  amplitude.  Only  the  very  initial  part  of  the  wave 
is  considered,  thus  the  pP  and  other  reflections  can  be 
neglected. 

The  recorded  amplitude  spectrum  T(f)  is  the  product  of  Ap(f) 
and  instrument  response.  It  has  a roughly  Gaussian  shape 
with  maximum  amplitude  A at  frequency  f , and  width  c. 

Ill  III  g 

Such  a spectrum  yields  on  Fourier  transform  a wave  packet 

with  amplitude  A = Amc//i  and  frequency  fm.  From  inspection 

of  recorded  spectra  it  was  concluded  that  a relationship 

c = f roughly  describes  the  spectra.  Since  the  distance 
m 

correction  in  m^  computation  compensates  for  the  Dp(A#h,E), 
we  get  for  m^ 

= log  F(fm)  + log  PQ  + ep  - 0.43  t%fm  + log|Wp| 

+ 2 log  fm  - log  /tt 

and  for  the  difference  m^  - Mg 

- Mg  = log  F(fm)  - 0.4  37it*  fm  + 2 Jog  fm  + C 

where  all  the  factors  which  do  not  depend  on  magnitude  are 
lumped  together  to  a constant  C.  They  are  either  constants 
or  random  quantities  with  some  expected  values. 

Using  the  same  source  spectral  models  and  f /m.  relationships 


at’  in  computation  of  the  predictions  for  the  LTM,  we  get  for 

m.  - M as  a function  of  m,  the  curves  which  are  shown  in 
d s d 27 

figure  VI. 2,  together  with  observations  by  Bungum  and  by 

Filson  and  Bungum28.  The  undefined  constant  C is  determined 

so  that  both  of  the  explosion  and  earthquake  curves  pass 

through  the  observed  points. 
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EXPLOSIONS  • 
EARTHQUAKES  on 


Fig.  VI . 2 . m.  - M values  of  explosions  in  Eastern  Kazakh  and 
earthquakes  in  Central  Asia.  Broken  lines  give  the  expected 
shapes  of  mb  - Ms  as  a function  of  m^,  levels  of  the  curves 
in  the  vertical  direction  are  determined  by  the  observations. 
mh  is  measured  only  at  NORSAR  for  events  with  square  symbols. 
Earthquakes  are  plotted  with  open,  explosions  with  filled 
symbols . 


The  essential  feature  of  the  two  curves  in  figure  VI. 2 is 
that  they  approach  each  other  as  magnitude  decreases.  The 
question  whether  they  will  converge  to  one  at  very  small 
magnitudes  was  treated  by  Douglas,  Hudson  and  Kembhavi  , whc 
showed  that  when  the  events  are  so  small  that  the  source 
time  function  appears  step- like  also  in  the  frequency  range 
of  the  short  period  sensor  response,  the  difference  mfa  - 
depends  on  which  of  the  following  two  factors  dominates: 

(a)  The  stronger  shear-wave  generation  of  the  double-couple 
earthquake  source  gives  for  the  quakes  larger  Mg  values  than 
for  the  explojions,  and  (b)  the  larger  depth  of  earthquakes 
decreases  the  earthquake  Mg  valuei  relative  to  explosions. 


The  contents  of  this  chapter  has  also  been  discussed  in 
reference  30. 
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VII  . 


SUMMARY  AND  CONCLUSIONS 


The  practical  aspects  of  deconvolution  of  a recorded  trace 
to  a record  of  earth  displacement  are  investigated  and  it  is 
concluded  that  the  urea  of  a unidirectional  wave  pulse  can  be 
determined  with  distortion  less  than  10  % if  the  pulse  lenght/ 
cut-off  period  ratio  is  less  than  0.35.  The  cut-off  refers 
here  to  the  filtering  away  of  the  lowest  frequency  components 
of  the  deconvolved  record,  containing  only  instrument  noise. 

E.g.  a pulse  of  2 seconds  in  lenght  requires  that  frequencies 
down  to  0.175  hz  be  included  in  synthesis  of  the  deconvolved 
trace,  using  FFT  methods. 

The  pulse  lenght  was  defined  as  the  time  interval  between 
the  pulse  onset  and  the  moment  when  the  displacement  returns 
back  to  its  pre-pulse  value.  The  pulse  width  was  defined 
as  the  pulse  area/height  ratio.  Pulse  width  is  usually  half 
of  its  lenght.  The  pulse  width  is  naturally  influenced  by 
the  amount  of  attenuation.  It  is  possible  to  correct  for 
some  attenuation.  A change  of  0.2  in  attenuation  parameter 
t*  produces  a change  of  20-30  % in  pulse  width.  For  explosions 
at  distances  less  than  30  degrees  from  NORSAR  often  no  distinct 
pulse  was  observed.  This  is  interpreted  as  caused  by  multipathing 

in  the  upper  mantle. 

The  common  spectral  model  consisting  of  a flat  low-f reouency 
model , a corner  frequency  fQ.  a high-frequency  constant  slope 
of  f'n  and  an  attenuation  factor  exp (-"ft  ) was  fitted  to  6 
explosion  and  earthquake  spectra,  grouped  regionally,  vhe 
t*-values  irom  all  studied  Asian  locations  to  NORSAR  were 
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observed  to  bs  less  then  0.4.  The  smallest  observed  attenuation 
(on  path  Kuriles-NORSAR)  has  t*=0.1  or  0.2,  depending  on 
assumptions  on  source  spectrum.  Values  (under  different 
assumptions  on  source  models)  for  the  ray  paths  from  Central 
Asia  to  NORSAR  range  from  0.15  to  0.25  in  t , corresponding 
to  a mean  Q equal  or  greater  than  2000.  Ray  paths  from  SW 
North  America  to  NORSAR  have  a larger  attenuation  than  any 

ray  path  under  Asia. 


The  simple  spectral  model  does  not  probably  quite  well 
describe  reality,  as  evidenced  by  an  unexpected  correlation 
among  the  derived  parameters,  and  in  some  cases  by  excessive 
corner  frequencies,  when  compared  to  pulse  widths  measured  by 
deconvolution.  A probable  deficiency  is  the  existence  of 
only  one  corner  frequency  In  the  model,  while  actually  several 
corners  may  appear,  separated  by  Intermediate  spectral  slooos 
Both  of  the  values  n*2  and  n=3  were  tried,  no  clear  preference 
for  either  one  was  deduced.  Corner  frequencies  could  be 
estimated  simultaneously  with  a regional  average  t , but  with 
a considerable  scatter.  Fitting  of  a spectral  model  to 
explosion  spectra  succeeded  only  if  a very  wide  frequency 
range  (up  to  6 hz)  could  be  used. 
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Pulse  widths  were  measured  for  the  same  events,  spectra  of 
which  had  been  studied.  Pulse  widths,  as  also  corner 
frequencies,  are  correlated  with  event  magnitude.  Explosions 
have  pulse  widths  about  half  of  those  of  earthquakes  with  same 
magnitude.  The  NTS  explosions  recorded  at  NORSAR  seem  to  have 
excessive  pulse  widths,  and  require  an  attenuation  correction 
larger  than  that  predicted  by  the  spectral  measurements  to  be 
reduced  to  the  same  pulse  width  as  the  E.  Kazakh  shots.  It 
is  important  to  bring  all  explosions  under  a same  statistical, 
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preferably  Gaussian,  population,  to  be  able  to  set  predictions 
and  statistical  limits  independently  of  region.  In  the 
following  table  observed  average  time  constants  of  seismic 
events  are  given 


Table 

VII. 1. 

Time 

constants . 

W=pulse  width  (s)  , 

duration,  f = 
o 

corner 

frequency 

(hz) 

earthquakes 

explosions 

% 

W 

D 

fo 

W fo 

4.0 

0.40 

0.8 

0.8 

5.0 

0.65 

1.2 

0.5 

0.25  1.3 

6.0 

1.1 

2.2 

0.3 

0.32  1.0 

The  pulse  width  seems  to  function  as  a discriminant.  It  is 
the  author's  opinion  that  after  the  m,  - M discriminant,  the 

D S 

pulse  width  is  best  in  separation  capability  and  applicability. 
It  was  tested  to  an  independent  sample  of  earthquakes  and 
explosions  recorded  at  LASA.  A benefit  of  the  deconvolution 
procedure  is  that  after  measuring  the  pulse  width,  one  is 
still  left  with  the  deconvolved  seismogram  which  frequently 
still  contains  useful  information  and  may  alert  the  observer 
in  case  of  multiple  events  etc.  Deconvolution  could  be  a 
standard  procedure  in  a monitoring  scheme  since  it,  as  broad- 
band seismograms  in  general,  displays  waveform  information 
easily . 


54 


APPENDIX 


In  the  following  table  the  first  84  events  on  a tape 
containing  LASA  beams  and  compiled  by  R.T.Lacoss  and 
J.R.Filson  (Lincoln  Laboratory,  M.I.T.)  are  listed. 

The  order  of  the  events  is  the  same  as  on  the  tape, 
the  epicentral  information  is  from  the  Bulletin  of  the 
International  Seismological  Centre,  or  in  a few  cases 
from  LASA,  in  these  cases  the  location  is  given  only  in 
full  decrees.  For  each  event  is  given  the  pulse  width 
measured  by  the  author  (P.W.)  and  the  low-frequncy  cut- 
off period  in  seconds  (C-0.).  When  the  cut-off/width 
ratio  is  less  than  4 and  according  to  the  discussion  in 
Chapter  II  one  may  assume  that  the  measured  pulse  width 
is  less  than  the  true  pulse  width,  the  word  "REDUCED"  is 

written . 
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FILE  1 


DATE  TIME  COORDINATES  DEPTH  MB  P.W.  C-O.  COMMENTS 


67-04-0-1 

12-23-35 

45. 7N 

151. 8E 

40 

5.9 

0.99 

:7.o) 

66-12-07 

17-17-42 

44. 3N 

151. 7E 

26 

5.8 

0.55 

7.0) 

67-10-07- 

•08-28-01 

49. 2N 

156.3E 

33 

5.3 

0.45 

>.o) 

67-11-01- 

■ 16-30-57 

48.  3N 

154.4E 

40 

5.5 

0.64 

:?.o) 

68-04-08 

10-10-24 

41. IN 

19.9  E 

4.2 

— 

:2.o) 

68-04-13 

01-39-04 

41  N 

22  E 

4.1 

— 

;2.o) 

68-04-11 

06-46-27 

42.  5N 

131. OE 

511 

5.0 

0.52 

Mo.) 

68-04-17- 

.iq_  14-00 

53. 9N 

160. IE 

33 

4.1 

— 

:2.o) 

68-01-27 

02-21-58 

40. 5N 

72. IE 

52  4.8 

— 

1.7) 

LASA  MB 

68-01-18 

21-57-03 

82. 3N 

119.5E 

48 

4.5 

0.62 

;i.7) 

REDUCED 

68-01-20 

10-42-08 

43. 2N 

147. 2E 

41 

5.2 

0.80 

M.o) 

68-01-20 

11-25-24 

43. 6E 

147. 2E 

33 

4.4 

0.36 

:2.o) 

FILE 

2 

67-11-50 

07-25-51 

41 . 5N 

20. BE 

29 

6.0 

0.45  1 

:i°.) 

67-12-25 

16-04-38 

48. 2N 

157. 5E 

26 

5.1 

0.45  ( 

:7.o> 

67-03-25 

18-62-17 

46. ON 

143. OE 

325 

4.8 



66-08- '9 

12-22-09 

39. 2N 

41.7  E 

26 

6.1 

1.94  I 

:7.o) 

68-04-  .7 

13_11_26 

36. 4N 

71.5  E 

113 

5.2 

0.91  1 

;7.o) 

68-04-18 

03-08-02 

41. 3N 

20. 3E 

33 

4.4 

0.43  1 

:2.5) 

68-04-21 

14-44-06 

54. 9N 

161. 5E 

28 

4.6 

0.42  1 

Cl. 7) 

68-04-17 

09-12-04 

35.2N 

3.7W 

16 

5.0 

0.54 

(7.0) 

68-01-29 

11-29-01 

42. 7N 

147. 6E 

33 

4.2 

0.40 

(1.7) 

68-01-29 

11-36-37 

43. IN 

146. 7E 

77 

4.2 

0.46 

(1.7) 

REDUCED 

68-01-29 

11-43-59 

43. 4N 

147. 3E 

33 

5.1 

1.08 

(3.0) 

REDUCED 

68-01-29 

12-07-08 

43. 2N 

147. 3E 

33 

5.0 

— — 

(1.7) 

FILE 

3 

67-05-25 

22-47-58 

45.  5N 

151. 4E 

41 

5.5 

2.0 

(7.0) 

66-05-20 

11-44-29 

55. ON 

165. 7E 

46 

5.2 

0.78 

(7.0) 

66-12-14 

1^-49-59 

45. 6N 

26. 4E 

158 

4.8 

0.27 

(4.0) 

67-01-05 

10-07-58 

39. 4N 

72. 9E 

11 

5.3 

0.64 

(2.5) 

REDUCED 

68-04-23 

02-43-45 

38. 7N 

70. 6E 

29 

4.8 

— 

(i.7) 

68-04-23 

22-30-27 

34. 6N 

8.9E 

33 

4.2 

0.41 

(1.7) 

66-10-27 

05-57-58 

73.  4N 

54. 8E 

0 

6.3 

clips 

68-07-01 

04-02-0 2 

47. 9N 

48.  OE 

33 

5.5 

0.28 

(4.0) 

68-01-30 

01-48-29 

43. 3N 

147. 7E 

33 

5.1 

0.58 

(i.5) 

REDUCED 

68-01-29 

17—14—06 

43. 4N 

147. 3E 

33 

4.5 

— 

(1.7) 

68-01-29 

16-42-50 

43. 5N 

147. 2E 

36 

5.7 

1.24 

(4.0) 

REDUCED 

^8-01-29 

16-04-12 

43.  IN 

146. 9E 

40 

4.6 

0.43 

(2.5)  LASA  MB 

i 


56 


FILE  4 


DATE 

67-01-11 

67-04-01 

67-01-20 

66-11-21 

66- 03-20 

67- 04-07 

67- 01-03 
66-05-07 

68- 01-30 
68-01-30 
68-02-06 
68-02-07 


TIME 


COORDINATES  DEPTH  MB  P.W.  C-0.  COMMENTS 


11-20- 

05-57- 

01-57- 

12- 19- 
05-49- 
^8-33- 
23-58- 

13- 08 
08-17 
01-30 
10-30 
01-22 


.46  34. IN 
09  46. 3N 
•23  48. ON 
.27  46. 7N 
■ 58  49. 7N 
• 31  37. 4N 
-21  48.  IN 
-16  37. 8N 
.32  36.4N 
i-l 3 43. 3N 

-34  54. 9N 
06  30. 9N 


45. 7E 
157. OE 
102. 9E 
152. HE 
78. OE 
36. 2E 
102. 9E 
27. 9E 
70. 7E 
146. 8E 
162. IE 
80. 3E 


34  5.6 
40  5.5 
33  6.1 
40  5.6 
0 6.0 
39  5.0 
33  5.4 
12  5.2 
205  5.2 
12  5.3 
33  4.2 
4.7 


0.77 
0.86 
1.17 
0. 56 
0.33 
0.96 
0.73 
1.28 
0.43 
0.73 
0.36 


g:Si 

7.0) 

4.0  ) 

7.0) 
4.0> 

4.0) 

7.0) 

2.5) 

4.0) 
(1.7) 
(1.7) 


67-04-01 

67-04-27 

67-05-27 

67-05-27 

67- 04-20 
67-02-26 
66-06-04 
66-10-29 

68- 02-07 
68-02-13 
68-03-17 
68-02-22 


14-00- 
23- 1 5 

01- 42- 
19-05- 

04- 07- 
03-57- 

05- 11- 

02- 39- 
22-22- 
07-42. 
17-46 
12-34 


34  45. 8N 
20  41. 7N 
•47  39. 9N 
■48  36.  IN 
>58  49. 7N 
■58  49. 8N 
■54  36. 3N 
-29  39. 2N 
-20  36. 7N 
-16  38. 4N 
-53  79. -yN 
-11  41. 5N 


FILE  5 

151.7E 
82.  3E 
77. 3E 
77. 8E 
78. IE 
78. IE 
70. 8E 
21. 2E 
26. 8E 
73. 2E 
140. OE 
20.  OE 


23  5.4 
33  5.0 
33  5.4 
35  5.4 
0 5.7 
0 6.0 
207  5.7 
20  5.7 
161  5.0 
33  4.3 
33  4.5 
33  4.1 


0.87 
0.59 
0.77 
0.93 
0.24 
0.29 
0.67 
0.49 
0.59 
0.47 
0. 56 


(4.0) 

(2.5) 

(1.7)REDUCED 

(4.0) 

(7.0) 


4.0) 

(7.0) 


(7.0) 

(7.0) 

(1.5)  REDUCED 
(1.7) REDUCED 

(1.5)  LASA  MB 


67-06- 

67-10- 

67-12- 

67-12- 

67- 01- 

68- 05- 
68-05- 
66-08 
68-02 
68-02 
68-01 
68-01 


07  18- 
05  15- 
14  IB- 
16  20- 
18  05- 
■ 11  12- 

• 13  02- 
•19  03- 
■22  12- 

• 22  09- 
-11  18- 
-08  09- 


^6-31 

55-03 

25-17 

53-58 

34-33 

12-41 

•46-36 

.53-OI 

22-50 

• 55-47 

-08-38 

-28-04 


47. 6N 
45. 4N 
54. 6N 
51. 2N 
56. 6N 
41. ON 
43.  5N 
50. 4N 
41 . ?N 
41. 8N 
46. 4N 
49. 7N 


FILE  6 

155. 4E 
150. 7E 
160.42 
157. 7E 
120. 8E 
49.  8E 
40. 3E 
77. 9E 
20. IE 
21. 5E 

153. 3E 
155. 9E 


29  5.2  0.58  (7.0) 

33  5.3  0.82  (2.5) 

33  5.5  0.73  (7.0) 

24  5.5  1.17  (7.0) 

11  6.1  1.  13  (7.0) 

16  5.0  0.43  (2.5) 

5 5.1  0.93  (7.0) 

0 4.8  0.21  (1.7) 

46  4.4  0.63  (2.5)  LASA  MB 
4.1  0.60  (l.5 )REDUCED 
50  4.7  0.23  (1.5) 

63  4.3  0.52  (1.7 )REDUCED 
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FILE  7 


DATE 

TIME  COORDINATES 

1 DEPTH  MB 

P. 

,w. 

C-0. 

COMMENTS 

67 

12- 

19 

i 

o- 

fO 

1 

o 

38 

N 

72 

E 

4.7 

0. 

.63 

Cl. 7) 

LASA  MB 

65 

11- 

21 

04-57-58 

49. 

8N 

78. 

IE 

0 

5.8 

0. 

.29 

C2.5) 

67 

10- 

30 

06-03-58 

49. 

8N 

78. 

IE 

0 

5.5 

0, 

,28 

C4.0) 

66 

05- 

07 

03-57-58 

49. 

7N 

77. 

9E 

0 

4.9 

0. 

,21 

1.7) 

65 

10- 

29 

21-00-00 

51. 

4N 

179. 

2E 

0 

6.1 

0. 

.37 

4.0) 

67 

04- 

07 

19-39-13 

47. 

ON 

146. 

0E 

296 

5.0 

0. 

,74 

4.0) 

68 

06- 

24 

18-2o-46 

40. 

4 A 

48. 

8E 

33 

4.6 

0. 

,38 

1.7 

68 

06- 

11 

03-05-58 

49. 

‘JN 

78. 

2E 

0 

5.3 

0. 

,31 

!2.5) 

68 

01- 

07 

21-38-52 

45. 

yK 

150. 

7E 

111 

4.5 

0. 

42 

:i.7) 

68- 

01- 

07 

18-45-48 

55 

N 

164 

E 

4.1 

0. 

48 

:i.7)  REJnfJCFin 

68- 

01- 

07 

07-42-04 

4-3. 

7N 

134. 

0E 

00 

0 

,-r> 

4.5 

mm  n 

. — _ 

;i.7) 

68- 

01- 

05 

02-51-33 

88 

N 

50 

E 

4.0 

0. 

50 

!i.7): 

REDUCED 

I 


REFERENCES 


1.  Randall,  M.J.,  "Spectral  peaks  and  earthquake  source 
dimension",  Jour.  Geophys.  Res.,  Vol.  78,  pp.  2609- 
2611  (May  1973) . 

2.  Haskell,  N.A.,  "Analytic  approximation  for  the  elastic 
radiation  from  a contained  underground  explosion",  Jour. 
Geophys . Res . , Vol.  12,  pp.  2583-2587  (May  1967). 

3.  Randall,  M.J.,  "Seismic  radiation  from  a sudden  phase 
transition",  Jour.  Geophys.  Res.,  Vol.  71,  pp.  5297- 
5302  (November  1966) . 

4.  Bogert,  R.P.,  M.J.  Healy  and  J.W.  Tukey,  "The  freauency 
analysis  of  time  series  for  echoes",  Time  Series  Analysis, 
edited  by  M.  Rosenblatt,  John  Wiley  & Sons,  Inc. (1963) . 

5.  Savage,  J.  , "Radiation  from  a realistic  model  of 
faulting",  Bull.  Seism.  Soc . Am.,  Vol.  56,  pp.  577-592 
(April  1966) . 

6.  Carpenter,  E.W.,  "Absorption  of  elastic  waves  - an 
operator  for  a constant  Q mechanism,  United  Kingdom 
Atomic  Energy  Authority,  AW RE  Report  D-43/66  (October 
1966) . 

7.  Frazier,  C.W.,  "Observations  of  pP  in  the  short-period 
phases  of  NTS  explosions  recorded  at  Norway",  Geophys . 

J.R.  astr.  Soc.,  Vol.  31,  pp.  99-109  (1972). 

8.  Aki,  K.,  "Scaling  law  of  seismic  spectrum",  Jour.  Geophys. 
Res. , Vol.  72,  pp.  1217-1231  (February  1967). 

9.  Randall,  M.J.,  "The  spectral  theory  of  seismic  sources", 
Bull.  Seism.  Soc.  Am.,  Vol.  63,  op.  1133-1144  (June  1973). 

10.  Berckhemer,  H.  and  K.H.  Jacob,  "Investigation  of  the 
dynamical  process  in  the  earthquake  foci  by  analyzing 
the  pulse  shape  of  body  waves".  Final  Sci.  Rept.,  1 June 
1964  to  31  December  1967,  Contract  AF  61  (052) -801, 

Inst,  of  Meteor.  & Geophys.,  Univ.  of  Frankfurt  (1968) , 


« 


i 


\ 

1 


59 


4 


Burdick,  L.J.  and  D.V.  Helmberger,  "Time  functions 
appropriate  for  deep  earthquakes",  Bull.  Seism.  Soc. 

Airu,  Vol . 64,  pp.  1419-1428  (OctobeF'l'S74)  . 

R,P*  a^d  s*s*  Alexander,  "Compressional  \ 'locity 
distribution  beneath  Scandinavia  and  Western  Russia", 
Geophys.  J.R.  astr.  Soc.,  Vol.  39,  pp.  587-602  (1974). 

Noponen,  I.,  "Compressional  wave  power  spectrum  from 
seismic  sources",  Sci.  Rept.  No.  1,  June  30,  1972  to 

i^j).Grant  OSR"72":^"7>  University  nf  Helsinki 

Haskell,  N. A. , "Total  energy  and  energy  spectral  density 
of  elastic  wave  radiation  from  propagating  faults".  Bull 
Seism.  Soc.  Am.,  Vol.  54,  pp.  1811-1841  (December  19*477’ 

Brune,  J.N.,  "Tectonic  stress  and  the  spectra  of  seismic 
shea r^ waves  from  earthquake*  , Jour.  Geophys . Res.. 

Vol.  75,  pp.  4997-5009  (Septembgr"TS76 ) . 

Savage,  J.C.,  "Relation  between  P-  and  S-wave  corner 
frequencies  in  the  seismic  spectrum",  Bull.  Seism.  Soc 
Anu,  Vol.  64,  pp.  1621-1627  (December  TT74)  . 

Savage,  J.r.,  "Relation  of  corner  frequency  to  fault 

o^ruuiy'^V^r'  Ge°phyS-  ReS->  Vol‘  77'  pp-  3788" 

Frazier,  C.W.  and  J.R.  ^ilson,  "A  direct  measurement 
of  the  Earth's  short  period  attenuation  along  a tele- 

3787™  (lg?^"'  — U-~-  Geophys-  Res->  Vol.  77,  pp . 3782- 

Seggern.  D.  von  and  R.  Blandford,  "Source  time  functions 
and  spectra  for  underground  nuclear  explosions",  Geophys. 
uS^-r • r>oc.  * Vol.  31,  pp . 83-97  (December  19  72"). 

Burridge,  R.  and  L.  Knopoff,  "Body  force  equivalents 
for  seismic  dislocations",  Bull.  Seism.  Soc.  Am,,  Vol. 

54,  pp.  1875-1888  (December  1964) . ~ 

Muller,  G. , Seismic  moment  and  long-pe.  od  radiation 
of  underground  nuclear  explosions" ,’ Bull . Seism.  Soc 
Am.,  Vol.  63,  pp.  847-857  (June  1973) . " 


60 


i 


21.  Randall,  M.J  . , "Low-frequency  spectra  in  seismic  waves 
from  explosions",  Geophys . J.R.  astr.  Soc.,  Vol.  32, 
pp.  387-388  (1973). 

22.  Douglas,  A.,  J.B.  Young  and  J.A.  Hudson,  "Complex  P- 
wave  seismograms  from  simple  earthquake  sources", 

Geophys.  J.R.  astr.  Soc.,  Vol.  37,  141-150  (1974)  . 

23.  Acton,  F.S.,  "The  classical  model",  Chapter  2,  Analysis 
of  Straight-Line  Data,  Dover  publ.,  Inc.  (1966). 

24.  Lacoss,  R.T.,  "A  large-population  LASA  discrimination 
experiment" , Technical  Note  1969-24,  Lincoln  lab., 

M.I.T.  (April  1965) . 

25.  Chinnery,  M.A.  and  R.G.  North,  "The  moment-Mo  relation- 
ship, and  the  frequency  of  large  earthquakes",  Semiannual 
Tech.  Summary  Rept.,  30  June  1974_,  Lincoln  Laboratory, 

M.  I.T.  (June  19  7 4)“. 

26.  Douglas,  A.,  J.A.  Hudson  and  V.K.  Kembhavi,  "The  analysis 
of  surface  wave  spectra  using  a reciprocity  theorem  for 
surface  waves",  Geophys.  J.R.  astr,  Soc.,  Vol.  23,  pp. 
207-223  (1971) . 

27.  Bungum,  H.  and  I.  Noponen,  "nu:M  at  NORSAR" , NORSAR 
Semiannual  Tech.  Rept.,  1 January  - 30  June  1974, 
ed.  by  H.  Bungum,  Kjeller  (September  19  74)  . 

28.  Filson,  J.  and  H.  Bungum,  "Initial  discrimination  results 
from  the  Norwegian  seismic  array",  Geophys.  J.R.  astr. 
Soc. , Vol.  31,  pp.  315-328  (December  '972) . 

29.  Douglas,  A.,  J.A.  Hudson  and  V.K.  Kembhavi,  "The  relative 
excitation  of  seismic  surface  body  waves  by  point 
sources",  Geophys.  J.R.  astr.  Soc.,  Vol.  23,  pp.  451- 
460  (1971). 

30.  Noponen,  I.,  "Event  identification  capabilites  of 
NORSAR" , Proceedings  from  NATO  A.S. I.  "Seminar  on 
exploitation  of  seismograph  networks,  Sandef jord, 

Norway,  1974  (to  be  published) . 


61 


